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Abstract: The origin of the ultra high energy cosmic ray (UHECR) show- 
ers has remained as a mystery among particle physicists and astrophysicists. 
In low scale gravity models, where the neutrino-nucleon cross section rises 
to typical hadronic values at energies above 10 20 eV, the neutrino becomes 
a candidate for the primary that initiates these showers. We calculate the 
neutrino-nucleon cross section at ultra high energies by assuming that it is 
dominated by the production of p-branes. We show, using a generalized 
Randall- Sundrum model, that the neutrino-nucleon cross-section at neutrino 
energies of 10 11 GeV is of the order of 100 mb, which is required for explain- 
ing UHECR events. Similar result also follows in other models such as the 
Lykken-Randall model. 
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Low scale gravity models [1, 2] lead to the possibility that the neutrino- 
nucleon cross section at ultra high energies may be several orders of magni- 
tude larger than what is predicted by the Standard Model. At center of mass 
energies ^/s less than the scale of quantum gravity M*, which we assume is 
of order 1 TeV, the cross section can be calculated using the perturbative 
Feynman rules [3]. However at higher energies, such that ^/s > M*, there 
does not exist any reliable procedure to calculate the cross section. In ap- 
plications to cosmic rays we are interested in cross section at center of mass 
energies of the order of 10 12 GeV 2 , where the perturbative approach is cer- 
tainly not applicable. Several authors [4, 5, 6, 7, 8, 9, 10, 11] have obtained 
various estimates by using different models for the cross section in this en- 
ergy regime. In Ref. [5] we assumed two models for the s dependence of the 
cross section, a ~ s, s 2 , and found that neutrino-nucleon cross sections are 
of the order of one to several hundred mb with the precise value dependent 
on the model used and on the choice of M*. The larger value of the cross 
section is obtained using the s 2 model. Such large values of the neutrino- 
nucleon cross section are very interesting since then the neutrino becomes a 
candidate [5, 6] for explaining the puzzling cosmic ray events which appear 
to violate the GZK bound [12, 13, 14, 15]. The idea that neutrino-nucleon 
cross section grows large at GZK energies is very old [16]. Large extra dimen- 
sion, low scale gravity models [1, 2] have provided new impetus to this idea 
[5, 6, 7, 17, 18, 19]. Even if the cross section is not large enough to explain 
the GZK events, cosmic ray observations can be used to test this hypothesis 
and to put stringent bounds on these models. 

There exist several proposals for calculating the ultra high energy scat- 
tering cross sections including the gravitational contribution. In Ref. [4, 9, 
20, 21] the authors proposed an eikonal model for the cross section, which 
may be applicable for very large s but for momentum transfer —t < M 2 . 
For larger values of the momentum transfer, Ref. [4, 22, 23, 24, 25] have 
proposed that black hole production will set in. This proposal has gener- 
ated great interest with several studies investigating the production rate of 
black holes at colliders [24, 26, 27, 28, 29, 30, 31] and cosmic ray collisions 
[10, 11, 32, 33, 34]. Using the eikonal model [9] and the black hole production 
rate [10], it is found that the neutrino-nucleon cross section is of the order 
of 0.1 mb at s pa 10 12 GeV 2 . This cross section is also large enough such 
that neutrinos could generate horizontal showers, which can easily be seen 
in future observatories. 
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Black hole production cross section is estimated by assuming that it is 
approximately equal to the geometric area of the black hole produced. A 
collision between partons % and j with center of mass energy is assumed 
to produce a black hole of mass Mbh ~ V§ and hence the parton level cross 
section a(ij — > BH) is equal to 7ri?| where Rs is the Schwarzschild radius of 
a black hole of mass Mbh- This geometric cross section has been criticized by 
Voloshin [35] on the grounds that a high energy collision with partons % and 
j will have a very large amplitude to radiate, and the amplitude to produce 
an isolated black hole will be exponentially damped. At our current stage 
of understanding of quantum gravity, the various models proposed must be 
considered as highly speculative. 

Recently Ref. [36] has proposed that besides black holes, high energy 
collisions will also produce branes. They argue that in certain cases the 
production cross section of p-branes of mass M p m V§ is much higher than 
that of black holes of the same mass. The precise value of the cross sec- 
tion depends on the dimensionality of the brane and the size of the extra 
dimensions. The largest cross section for a p-brane is obtained when the 
brane is completely wrapped on the small-size extra dimensions where it is 
assumed that there exist m extra dimensions compactified on length scale of 
order L<M~ l and the remaining n — m extra dimensions are compactified 
on length scale of order L' » M~ l . The ratio E(s; n,m,p < m) of the cross 
section of the p-brane to black hole production in this case is given by [36] 

r ~ 2p / \9 w-l 

/ L \ "-P+1 7(71, p) ( S \ 

E( S ; W < m) «(-) ^(-) (1) 
where n is the total number of extra dimensions, = M" 1 , 

ST (^±p) 



7(n,p) = 



n+l 



(2) 



and w — 1/[1 — p/(n+ 1)]. The branes, once produced, will decay producing 
gravitons and matter fields. 

We calculate the neutrino-nucleon total cross section a up at ultra high 
energies assuming that the cross section is dominated by brane production. 
We perform this calculation within the framework of a generalized Randall- 
Sundrum (RS) model [2] and the model proposed by Lykken and Randall in 
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Ref. [37] such that these models have an arbitrary number of compact extra 
dimensions instead of the single compact dimension assumed in Ref. [2, 37]. 
Besides the single compact dimension which involves the warp factor, the 
remaining m compact dimensions are assumed to be factorizable. These m 
extra dimensions are compactified on a m-torus of radius R. The generalized 
RS model, therefore, consists of two 3-branes with opposite tension situated 
at orbifold points y = 0, y c in a (4+m+l)-dimensional space-time. The metric 
can be written as 

ds 2 D+1 = e~ 2a{y) r 1 ^dx^dx y + R 2 5 ab d9 a d9 b + dy 2 (3) 

where /i, v = 0, ..,3, a,b = l,...,m, 6 a are the compactied dimension with 
range [0,2n] and y is the warped dimension with range [-y c , y c ]. The Einstein 
equations can be written as 

Gab = 4M ™+3 ( T *b + T ab) ( 4 ) 
The components of the bulk energy-momentum tensor are assumed to be 

T^ = ^A , 7f>* = Ao, T^ b = A e 5 b a 
The components of the brane energy-momentum tensor are assumed to be 



'm i 



where i=l,2, and all other non-diagonal components are taken to be zero. In- 
homogeneities in the energy-momentum tensors and cosmological constants 
can be due to different contributions to casimir energies in the different direc- 
tion of space-time [38, 39] . In the absence of a four dimensional cosmological 
constant, the (//, z/)-component of the Einstein's equation gives 

-3o"(y) + Qa-(y) = - £ ^S(y - Vi ) (5) 

i 

The (yy)-component of the Einstein's equation is given by 



4 



The (9, ^-component of the Einstein's equation is given by 



-4o"(y) + 10.'%) = - £ ^6(y - Vi ) (7) 

i 

From Eq. (5,6,7), using the orbifold symmetry along y-direction, we find 



a(y) = V~ 24M°™+3 ^~ Vc) (8) 

and A e = |A , i$ = |i>q- This solution is valid only if brane tensions and 
cosmological constants are related as follows 

A = -24M m+3 k 2 , vl = -v 2 = 24M m+3 k (9) 
The four-dimensional reduced Planck mass M v \ can be written as 
M 2 = M (m + 3) exp (2%c) (2^T fl _ exp( _ 2 % c ))] 



(10) 



where M = ( 327r )i/(™+3) ■ If we demand ky c ~ 35, M* = lTeV and R < 
(TeV)~ l we find that k < M*. We can take y c of order fm without violat- 
ing the existing observational constraints from astrophysics and accelerator 
physics and then we find k is of order GeV. We assume that some unknown 
mechanism stabilizes the radius of extra dimensions. 

We next determine the spectrum of the KK-modes of this theory in the 
low energy limit. Decomposition of KK states can be carried out by consid- 
ering the following graviton perturbations: 



ds 2 = e~ 2a \^ v + h^(x, y, 9^} + dy 2 + R 2 5 ij d9 i d9 J (11) 

where R is the radius of compact extra dimensions and a = k(\y\ — y c ). Here 
we are not considering fluctuations around extra dimensions. We expand the 
linear fluctuations around the flat 4-dimensional metric in a complete set of 
radial wavefunctions and fourier modes: 



K„{x, v a) = T, h ^Hx)<P(ndyy w (12) 
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Working in a gauge if^ 'd a h^ = t] al3 h^ = and d 2 hfy l \x) = m 2 Ji^J^{x) 
the differential equation for radial part is given by 



d y (e-^dy^M) = [ml - W2 e- 2 °) My) (13) 



with normalization condition (2irR) m J^_ c dye 2a <fi* n 4> n i = 8 nn >. After chang- 
ing variables to z n = ^e a and f n j = e~ 2<T ni ;, Eq.(13) can be written as (for 
y ^ 0, y c ) 



„2 d2 fn,l + z dfn£ + 



k 2 R 2 



fn,l = (14) 



In the low energy limit we can set 1=0 since these modes will have energy 
larger than 1/R. Then in this limit the theory reduces to RS model except 
for the fact that the size of the warped extra dimension is taken to be of order 
1 fm which is about three orders of magnitude larger than the fundamental 
length scale of (ITeV)" 1 . The masses of these KK-modes are quantized in 
mass scale of the order of GeV and their couplings to the standard model 
fields which lives on the brane at y = y c are the same as in RS case. As 
there are a large number of modes present, we can ignore the few low lying 
modes and the sum over the remaining modes can be done by changing the 
summation to an integration as in the ADD case. 

We also consider another model proposed by Lykken and Randall [37] 
in which the size of the extra dimension can be infinite. The model was 
originally proposed with only one extra dimension. The model involves three 
branes, one located at y — 0, second at y = oo and the third, which is the 
physical brane, located at some position y = y such that < y^ < oo. We 
generalize this model by introducing m additional compactified dimensions 
in exact analogy to the above construction in the case of RS model. We find 
that the low energy predictions of the model remain unchanged due to the 
addition of m compact extra dimensions. Hence we find that we can extend 
both the RS and the Lykken-Randall model such that they contain m extra 
dimensions compactified on a length scale of the order of inverse TeV and 
one large extra dimension of length scale larger than (or order of) 1 fm. We 
can now construct p-brane solutions in these models such that the brane is 
completely wrapped around the m extra compact dimensions. 

The parton level brane production cross section o"b ra ne for p-brane of mass 
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M p = V§ is given by, 

o-branc(s) = £(s;n,m,p < to)(7 B h(s) (15) 



where <5"bh is the production cross section for black hole of mass Mbh = 
The black hole production cross section, <3"bh, is given by [23, 24] 




where rs is the Schwarzchild radius of a 4 + n dimensional black hole, 



The black hole and the brane production processes are expected to give 
dominant contributions when s » M%. 

The total cross section can be computed using 



where is the parton distribution function corresponding to the i parton, 
x m i n = Ml j s and Q is the typical momentum scale of the collision which is 
taken to be the brane mass. We use the CTEQ parton distributions [40] for 
our calculation. Since the parton distributions are not known beyond Q — 10 
TeV, we set Q equal to this value if the brane mass exceeds 10 TeV. The 
minimum value of x is obtained by assuming that branes are produced with 
masses M p greater than M*. 

In fig. 1 we plot the neutrino-nucleon cross section for some representative 
choice of parameters, p and n. The ratio of length scales L/L* has been set 
equal to 0.25. We find that for n = 7 and p = 6 the cross section rises 
to roughly 50 mb at neutrino energy of 10 11 GeV. As p becomes smaller 
the cross section is much smaller. Similarly the cross section is considerably 
reduced as we lower n. The cross section is also found to be very sensitive to 
the precise value of L/L*. The dependence of total cross section on L/L* is 
shown in fig. 2 for some representative choice of parameters, p and n. The 
neutrino energy has been taken to be 10 11 GeV for this plot. 

It is clear from fig. 1 and fig. 2 that there exists a large range of parameter 
space where the neutrino-nucleon cross section is of the order of 100 mb or 




(17) 




(18) 
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Figure 1: The neutririo-nucleon cross section a vp including the p-brane pro- 
duction for several representative values of p and the number of extra di- 
mensions n. The ratio of the length scales L/L* of the small to large extra 
dimensions has been set equal to 0.25 and the scale of quantum gravity 
M* = 1 TeV. The highest energy HERA data point and the cross section 
obtained in Standard Model (SM) are also shown for comparison. 
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above for the primary neutrino energy of 10 11 GeV. These are comparable 
to the nucleon-nucleon cross section at these energies and is typically the 
cross section values required in order that the neutrino can considered as a 
candidate for generating the observed air showers above the GZK bound. We 
point out that in the present case the neutrino delivers its entire energy in a 
single collision. This is in contrast to our earlier work [41], using a t-channel 
exchange in which the neutrino loses only a fraction of its energy in a single 
collision. In that case, the neutrino collides with several air nuclei along its 
path, and the resulting air shower tends to be longer in developing compared 
to an equally-energetic proton shower. Nevertheless, the characteristics of 
air showers do not rule out this model, given fluctuations and uncertainties 
in energy [41]. In the present case with high inelasticity, the position of 
the shower maximum will be determined dominantly by the value of the 
neutrino-air cross section. Since there exists a large range of parameter 
space where the neutrino-nucleon cross section exceeds 100 mb, the shower 
maximum could exist at the same or higher altitude than that produced by 
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Figure 2: The L/L* dependence of the neutrino-nucleon cross section a up 
for several representative value p and the number of extra dimensions n. he 
scale of quantum gravity M* = 1 TeV and the incident neutrino energy is 
taken to be 10 11 GeV. 
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a proton or even an iron nucleus of equivalent energy. Hence the position 
of the shower maximum may be indistinguishable from a hadron primary of 
equivalent energy. 

We point out that our cross section does not violate any bounds imposed 
by unitarity. The models and parameter ranges considered in this paper 
are not ruled out by the low energy data. The cross section estimate is 
based on the calculation of the area of the classical brane solution, which 
essentially can be regarded as an estimate of the size of the region where 
the gravitational interaction between the two particles is strong. Therefore 
within these models, the total cross section has to be atleast as large as the 
estimate obtained by assuming that it is dominated by p-brane production. 
Furthermore as argued in Ref. [41] it is not possible to determine whether a 
cross section violates unitarity simply by looking at its rate of growth with 
energy. For example, the large o v - v (> lOOmb) obtained in Ref. [41] by 
assuming that it rises as s 2 with energy, is well within the bounds obtained 
by Goldberg and Weiler [42] by a model independent application of dispersion 
relations and unitarity. The growth of p-brane production cross section with 
energy is relatively mild in any case, with the total cross section rising at most 
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as E^ 6 at large energy for the parameter ranges considered in this paper. 

After formation the p-brane will decay. In analogy with black hole evap- 
oration, we expect that branes will decay into gravitons and matter fields. 
The emitted gravitons will not be detected. Hence only a fraction of the 
particles that are produced will be responsible for generating the shower and 
the primary energy may be significantly underestimated. This implies that 
for a shower of energy E we need a primary neutrino of energy E' which 
is somewhat larger than E. The precise amount of energy which goes into 
gravitons is model dependent. 

In conclusion we have shown that in low scale gravity models the neutrino- 
nucleon cross section is of the order of or larger than 100 mb for a large 
range of parameter space allowed by current experimental constraints. This 
makes the neutrino a candidate for explaining the cosmic ray events which 
appear to violate the GZK bound. Future cosmic ray observatories will be 
able to confirm or rule out this hypothesis. In particular, the ability of 
neutrino events to point back to their sources, which can be established by 
well-defined statistical correlations [43], has the potential to establish or to 
rule out neutrinos as primaries in the GZK-violation mystery. A very clear 
discussion on the possible origin of the highest energy cosmic rays is given in 
Ref. [44]. 
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